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Abstract 

Instanton matrix modeis (IMM) for two dimensionai string theories are obtained from 
the matrix quantum mechanics (MQM) of the T-dual theory. In this paper we study the 
connection between the IMM and MQM, which amounts to understand T-duality from 
the viewpoint of matrix models. We show that type OA and type OB matrix models 
perturbed by purely closed string momentum modes (or purely winding modes) have the 
integrable structure of Toda hierarchies, extending the well known results for c = 1 string. 
In particular, we show that type OA(OB) MQM perturbed by momentum modes has the 
same integrable structure as type OB (OA) MQM perturbed by winding modes, which is 
a nontrivial check of the T-duality between the matrix models. The MQM deformed by 
NS-NS winding modes are used to study type string in 2D black holes. We also find 
an intriguing connection between the IMM and the MQM via tachyon condensation. The 
array of alternating D-instantons and anti-D-instantons separated at the critical distance 
plays a key role in this picture. We discuss its implications on sD-branes in two dimensional 
string theories. 
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1. Introduction 

Recently the c = 1 matrix quantum mechanics(for reviews see [0,^,0]) has received a 
lot of attention because of its new interpretation as the decoupled world volume theory 
of unstable D0-branes[f|,|5],|6[]. The matrix models dual to type string theories were also 
proposed in |7|||. For other recent developments, see P, p^ , |lllP^Jl^ , |l^ , p!5l , |l^[l7| . The 
type OB matrix quantum mechanics (MQM) describes open string tachyons living on the 
unstable DO-branes, whereas the type OA MQM describes tachyonic open strings stretched 
between stable DO- and anti-DO-branes. Upon compactification on Euclidean time, these 
two matrix models are conjectured to be T-dual to each other. The exact agreement in free 
energy was found in [Q. However, unlike c = 1 matrix model which can be derived from 
discretizing the Polakov action on the string world sheet, such a derivation is not known 
for type matrix models. From the matrix model point of view, the T-duality between 
type OA and OB strings seems rather mysterious. To understand this T-duality is one of 
the motivations of this paper. 

To start, let us consider type OA MQM with Euclidean time compactified on a circle 
of radius a' / R. This is supposed to be T-dual to type OB string theory on a circle of radius 
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R. By decomposing the fields in type OA MQM into their Fourier modes along the thermal 
circle, we get a zero-dimensional matrix model, of the form 

/ dUdU J^^ n ^t e -/3Tr[(Xt„-t„X+27rnfit„)(Xt„-t„X+27rni?t„) t -a 2 t„tt] ^ -Q 

J n 

where U = e %x t R , U = e xX l R are the holonomies of the Wilson lines in type OA MQM, 
t n , £T are the Fourier modes of the complex tachyons. This is an instanton matrix model 
(IMM) , since the Wilson lines on the DO-branes in type OA theory are mapped to collective 
coordinates X, X for D-instantons in type OB theory, and the tachyons on the DO-branes are 
mapped to tachyonic open strings stretched between D-instantons and anti-D-instantons. 
Similar D-instanton matrix models have been studied in the context of ten dimensional 
type IIB string theory [[TBljnfl . Now having two different matrix model duals of type OB 
theory, we want to understand how they are related to each other. Since one of them 
involves unstable DO-branes, and the other involves D- and anti-D-instantons, it is natural 
to suspect that tachyon condensation plays the key role of connecting the two theories. If 
we can show the equivalence between the IMM and MQM, we would prove the T-duality 
between type OA and OB matrix models. 

The first thing we want to understand is the identification between operators in IMM 
and in MQM, in particular the ones that correspond to closed string momentum and 
winding modes. This was well understood in the context of c = 1 string theory. The 
momentum modes in MQM are represented as asymptotic perturbations of the fermi sea, 



whereas the winding modes are identified with the holonomy of the gauge fields p0| , |21 
c = 1 string perturbed purely by momentum modes []22"|j23|1 , or by winding modes has 
the integrable structure of Toda lattice hierarchy, with the integrable flow generated by the 
corresponding closed string perturbations. In particular, the perturbed grand canonical 
partition function is shown to be the r-function of the corresponding integrable hierarchy. 
The integrable structure appearing in c = 1 string theory is subject to constraints, known 



as the string equation ||24|| . These constraints can be equivalently thought of as imposing 
an initial condition on the flow of r-functions, which is the unperturbed partition function. 
Since the unperturbed c = 1 MQM on radius R and a'/R have the same free energy, it 
follows from the integrable structures that the grand partition function of two theories 
perturbed respectively by momentum modes and winding modes also agree. 



We shall generalize this approach to type string theories. In type OB MQM, the 
symmetric and antisymmetric perturbations of the fermi sea decouple.0 They generate two 
independent Toda flows, subject to different string equations. Consequently the perturbed 
partition function is the product of two r-functions, associated with the symmetric and 
antisymmetric perturbations respectively. On the type OB IMM side, one can integrate 
out the tachyons and get a unitary matrix model only in terms of the "holonomies" , or the 
collective coordinates of the instantons. This is T-dual to the "twisted partition" discussed 
in pO| , pTf in the case of c = 1 string. As we will see, the perturbed grand canonical partition 
function of the IMM indeed has a similar structure, provided a nontrivial identification of 
the perturbation parameters with those of MQM. We will use it to identify the operators 
in IMM that are dual to NS-NS and R-R scalars in spacetime. 

Type OA MQM and IMM perturbed by momentum modes also have the structure of 
Toda hierarchy. This is very similar to the case of c = 1 string, except that the string 
equations for type OA MQM and IMM are nonperturbatively well defined. The combination 
of these shows that both type OA and type OB theories perturbed purely by momentum 
modes, or purely by winding modes, are integrable. This also directly verifies that the 
IMM and MQM are equivalent at least when only momentum modes or winding modes 
are present. The perturbations involving both momentum and winding modes are more 
complicated, since in that case we would lose integrability. 

An alternative attempt to connect the IMM with MQM, which is more direct and 
intuitive, is via tachyon condensation. Consider turning on an open string tachyon profile 
on the DO-brane world volume T(X) ~ A cos(X/\ / 2a') 1 where X is the Euclidean time. 
As well known [p6| , |27|| this corresponds to a marginal deformation in the worldsheet CFT. 
With sufficiently large A, it takes a DO-brane into an array of alternating D- and anti- 
D-instantons separated at the critical distance. It is natural to expect that, the MQM 
expanded near this tachyon profile should be the same as the IMM expanded near the 
configuration of an array of D-instantons. 

On the MQM side, this periodic tachyon profile effectively discretizes the Euclidean 
time circle to a periodic lattice of spacing ny/^a'. It is well known]]]] that the MQM on 
a time lattice of spacing e with e less than a critical distance a is exactly equivalent to 
the continuum theory, provided proper redefinition of the parameters. In our case a = 



1 As remarked in |]2ql , this doesn't mean that the NS-NS and R-R scalars decouple, because 
of the nonlinearity of bosonization. 
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its/lot! is the same critical distance at which the tachyonic open string stretched between 
D- and anti-D-instantons becomes massless. The array of D-instantons is classically a 
stationary configuration in the IMM. If we integrate out the (complex) tachyons connecting 
the instantons, we expect an instability that drives the D- and anti-D-instantons toward 
each other. On the other hand, the collective coordinates of the instantons are eigenvalues 
of Hermitian matrices in the IMM. They effectively repel each other and fill up a "sea" of 
D-instantons. We have a large number of D-instantons distributed in a periodic effective 
potential V(X) along the thermal circle, and the D-instanton array corresponds to a critical 
point where the instantons are sitting at the top of the potential. There is a phase transition 
when the "instanton sea" merges the top of the potential, and this is the critical point that 
defines the double scaling limit of the IMM. This is analogous to the case of c = 1 matrix 
model, where the double scaling limit is defined as the limit that the fermi level approaches 
the top of the tachyon potential. By relating the collective coordinates of the array of D- 
instantons to the open string tachyons in the discretized MQM at the sites on the time 
lattice, we will show that the two matrix models become the same in the limit of critical 
distance. 

Fig. 1: Schematic picture of the array of D-instantons condensing along the 
Euclidean time circle into the "instanton sea" . 

It was proposed in |28| (see also [p9| , |30|| ) that sD-branes (A = 1/2 s-brane) are described 



by an array of D- and anti-D-instantons along the Euclidean time, separated at the critical 
distance (see also [PD[ ). In the sense described above, the IMM can be thought of as 



the world volume theory of sD-branes, and sD-branes play the same role in the IMM as 
DO-branes in the MQM. We compute the closed string fields sourced by sD-branes from 
IMM, and find exact agreement with calculation using ZZ boundary states. In fact, we can 
reproduce the (1, 1) ZZ boundary state in Liouville |31]] and super-Liouville theory | f32| , |33 | 
from c = 1 IMM and type OB IMM in a very simple manner. 
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The c = 1 matrix model deformed by the lowest winding modes is conjectured to 



describe bosonic string in a 2D Euclidean black hole [21]. The conjecture is extended to 



type theory in [34]. The exact free energy of the deformed type matrix models can in 



principal be obtained by solving the Hirota differential equations. Using the technique of 
|21||, we can solve for the genus and genus 1 terms in the perturbative expansion of the 



free energy. They are of the form[21,35 



T= -2Tv(R-R H )M + (s 1 + l)lnM + --- (1.2) 

where M, depending on the coefficient A of the winding mode deformation, is interpreted 
as the mass of the black hole, and Rjj is the asymptotic radius in Euclidean time. In the 
T-dual Af = 2 Liouville theory ]36|,[37| , A is essentially the coefficient of the Liouville term. 
We find from the genus piece of the free energy that M oc A 4 , which agrees with the 
general behavior expected from M = 2 Liouville theory. The density of states for the black 
hole has Hagedorn growth 

p{M) ~ M Sl e^ M (1.3) 

We find that s\ = —17/12 for type OB string and si = —13/12 for type OA string (un- 
charged black hole). 

This paper is organized as follows. In section 2 we derive the IMM for c = 1 and type 
strings. In section 3, we review the integrable structure in the S-matrix of c = 1 MQM, and 
extend them to the case of type OA and type OB MQM. We show that the same integrable 
structures appear in the corresponding IMMs, and that they are subject to the same string 
equations by computing the unperturbed free energies. Section 4 studies the connection 
between IMM and MQM from the viewpoint of tachyon condensation. In section 5, we 
compute the closed string fields sourced by sD-branes from IMM and compare them to the 
results obtained from ZZ boundary states. The integrable structures of type theories are 
applied to computing the free energy of the 2D Euclidean black hole in section 6. 



2. The Instanton Matrix Model 

In this section we derive the instanton matrix model of c = 1 string and type string 
theories from the (gauged) matrix quantum mechanics of the T-dual theory. We will set 
a' = 1 in the case of c = 1 string, and keep a' explicitly in most of the discussions on type 
theories, unless otherwise indicated. 
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2.1. The IMM for c = 1 string 

Consider the gauged c = 1 MQM with Euclidean time x compactified on a circle of 
radius R' = 1/R. The matrix model action is 

i(d x $ + z[A$]) 2 -^ 

where the tachyon and the gauged field A(x) are Hermitian N x N matrices. Let us 
fix the gauge d x A = 0, which sets A to its zero mode A^ and introduces the Fadeev-Popov 
determinant 



dxTr 



(2.1) 



db dcexp(Trbd x D x c) = J [ 



/ sin[(ai — a,j)/2R] \ 



-V (a.-a^^R J 



[2.2) 



where a^'s are the eigenvalues of A^°\ Let us decompose into its momentum modes 
$ < - n - ) along the thermal circle. The action is now written as 



S= J[J2 Tt \ (^ (n) + l Ai0) i $(n) ]) {nR^ - [A<°>, $ 



(-»)■ 



_ I<&W<b(-«) 
2 



(2.3) 

In the T-dual theory X = 27rA(°) is the collective coordinate for D-instantons, $* n * ) are 
the open string tachyons stretched between D-instantons with (relative) winding number 
n. T-duality acts on the inverse "Planck constant" (3 as (3 — > PR. We end up with an 
instanton matrix model action 

£ 

2tt 



-i(27r) 2 $ (n) $ ( " n) 



(2.4) 



The path integral measure is modified by the FP determinant ( |2.2|) . If we diagonalize X 
in terms of its eigenvalues Xi, the measure factor involving X becomes the measure for 
unitary matrices 

II sin2 (^r) ( 2 - 5 ) 

In other words, the natural variable to be integrated over is the "holonomy" U = e lX l R . 
As well known [ p8| , |39| , this can be thought of as the contribution from D-instantons at X{ 
together with all of their images at Xi + 2 r nnR. 

In the infinite radius limit R — > oo, we can drop all the modes $( n - ) with nonzero 
winding since they become infinitely "massive" . The action becomes simply 



S = £tt([X,^} 2 - 4n 2 (^y 
Arc 



(2.6) 



An alternative instanton matrix model for c = 1 string, the Liouville matrix model, 



was studied in [4C]. It would be nice to understand its relation to the IMM we are proposing. 
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2.2. The IMM for Type String Theories 



Let us start with the MQM of type OA theory in two dimensions, which is the de- 
coupled world volume theory of (stable) DO- and anti-DO-branes. In the background with 
no net DO-brane charges, the matrix model has U(N) x U(N) gauge symmetry. This 
is the case we will be concerned with. We have the U(N) x U(N) gauge field Aq and 
bifundamental tachyon (f>, 



A 0' 
I 




The Lagrangian is 



L = Tr 



(D ty (D Q t) + ^th 



(2.7) 



(2- 



where Dot = dot + At — t A. Again, we want to compactify the Euclidean time on a circle of 
radius R' = a' / R, and rewrite the matrix quantum mechanics in terms of a matrix integral 
over the Fourier modes of A and t. The DO- and anti-DO-branes becomes D- and anti-D- 
instantons; the zero modes of A and A (Wilson lines) become the collective coordinates 
of the D- and anti-D-instantons; the Fourier modes of t and v are open strings stretched 
between D- and anti-D-instantons, which can also wind around the circle. We expect this 
model to be equivalent to type OB MQM at radius R in the double scaling limit. 

Let us fix the gauge OqAq = 0, which sets A and A to their zero modes A^ = X/2ira' 
and A^ = X/2na'. As before the gauge fixing introduces the FP determinant 



-j-j- f sm[( Xl - Xj )/2R] \ 2 r S m[(x l -x J )/2R} \ 
11 V (x l -x J )/2R ) V (£i-%)/2i2 ) 



(2.9) 



where X{ and x% are the eigenvalues of X and X respectively. In terms of the Fourier 
modes t^ of t(x), the action is 







7r(2a') 3 / 2 



Tr [(2TvnRt {n) + Xt^ n) - t^X \ (^imRt^ + t (n)t X - A¥ n)t 



(2.10) 

where we exhibited the (type OB) inverse "Planck constant" /3, a = tx\J2ol' is the critical 
distance. There was a factor 27r a' / R coming from the integral over Euclidean time, but 
under T-duality (3 — ► f3R/\ / 2a' [§, so the factors of R cancel out. If we further gauge fix 



X to the diagonal form, the usual measure factor A(x) 2 A(x) 2 is converted to the measure 
for unitary matrices 

TTsin 2 ( — — — ) sin 2 ( — — — ) (2.11) 
11 v 2R 1 K 2R 1 K ! 

Therefore the natural variables to be integrated over are the "holonomies" U = 
e 2inX/R^ jj _ e 2-KiX/R^ rp^g (iiff erence disappears in the infinite radius limit R — > oo. 

In this limit all the winding modes are very massive and we can drop them, so the action 
simplifies to 

\xt-tX)(tfX -Xtf) -a 2 tt f (2.12) 



' TV 



7T 



(2a') 3/ 



Even though in the decompactification limit the integral over X, X is simply a Gaus- 
sian, we shouldn't integrate out them directly. One reason is that, the determinant coming 
from integration over X, X would completely cancel the F-P determinant coming from di- 
agonalizing t. We would then naively conclude that the eigenvalues of t do not repel each 
other, and it is unclear how to define the double scaling limit. What we have done wrong 
is reminiscent to the case of gauged MQM, where integrating out the gauge fields naively 
cancels the Vandermonde determinants of the eigenvalues at every time. This is incor- 
rect, since if we carefully discretize the Euclidean time, the gauge fields appear as "link 
variables", while the tachyon fields are associated with the "sites". At the end, the F-P 
determinants of the eigenvalues at the initial and end points of the time evolution are not 
completely cancelled out. Their effect is nothing but to antisymmetrize the wave functions. 
Similarly in the IMM before integrating out X and X, we should either restrict to a finite 
interval of time, or compactify Euclidean time on a finite circle. 

There is another way to understand this subtlety. An unstable DO-brane with A = 1/2 
rolling tachyon profile is an array of alternating D-instantons and anti-D-instantons along 
the thermal circle, separated at the critical distance a = twJIol' |2^1 - Suppose 2ttR = 2ma 
for some integer m. Naturally we expect that the U(N) MQM to be equivalent to the 
IMM with gauge group U(mN) x U(mN). Therefore in the limit R — > oo, we must 
correspondingly take iV — > oo. We will come back to this point in section 4. 

Clearly we can interchange the role of type OA and OB MQM in the above discussion. 
We will then obtain an IMM of type OA theory. Formally this model has the same action 
as the IMM for c = 1 string. Presumably the double scaling limits are defined differently 
in these two models. 
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3. Integrable Structures in Two Dimensional String Theories 

It is well known that c = 1 string deformed by closed string momentum modes (or 
winding modes) have the integrable structure of Toda chain hierarchy p^ , p4| , |23| , ^T| . The 



operators corresponding to the momentum modes generate Toda flows, and the perturbed 
grand canonical partition function is identified with a r-function. In this section we will 
show that similar integrable structures appear in type OA and OB string theories. Sub- 
section 3.1 is a review of some results on the exact S- matrix and integrable structures of 
"theory I" (c = 1 matrix model), essentially following |2^,^2[. There is nothing new in this 



first subsection, but it sets up the conventions that will be used in the rest of the section. 
In subsection 3.2, we study the S-matrix and integrable structures of "theory II" (type OB 
MQM). Most of these results are already known, but we will derive the string equations 
for theory II explicitly. In subsection 3.3, we show that type OA MQM also has the Toda 
integrable structure. This integrable structure is, in some sense, more natural than that of 
"theory I" since type OA theory is non-perturbatively well-defined. Subsection 3.4 studies 
the perturbed partition function of type OB IMM. We will indeed recover the structure 
similar to that of type OB MQM, and find a dictionary translating closed string modes 
in spacetime to operators in the IMM. In subsection 3.5, we compute the free energy of 
IMM by explicitly performing the matrix integral. The case for c = 1 IMM was essentially 
done in [[20,21]. In the case of type OB IMM, we find that the free energy factorizes into 



two pieces, which precisely agree with the contribution to type OB MQM free energy from 
symmetric and antisymmetric fluctuations of the fermi sea. By matching the unperturbed 
free energies, we conclude that the string equations for the integrable structure of IMM 
are the same as those for MQM. 

3.1. "Theory I" 

It is convenient to define light cone variables 

In "theory I" (c = 1 string) these variables are restricted to the region x± > 0, and for now 
we will not worry about nonperturbative effects. The Hamiltonian of the free fermions is 

Ho = — {x + X- + X-X + ) (3.2) 
9 



Since x+,x_ are conjugate variables, we can write the wave function of a state in either 
x+ or X- representation. The two wave functions i/j+(x+) and i/j-(x-) are related by a 
Fourier transform 

V>-0_) = (Sij) + )(x + ) 

r™ (3.3) 
= / dx + K(x-, x + )^fj + (x + ) 
Jo 

where the integration kernel K is 

[2 

K(x_,x + ) = W-cos(x_x + ) (3-4) 

We could have also chosen K(z-,z + ) = a/2/7t sin(z_z_|_) instead. This ambiguity reflects 
the fact that perturbatively the two sides of the fermi sea decouple. The energy eigenstates 
have wave functions 

, TP , n 1 HE— 77 

^(x±) = -=x ± 2 (3.5) 

V Z7T 



Note that we have chosen a convenient basis for i/j e and tpE, but they are not related by 
(p.3|) . In fact, it is straightforward to show that 



^Vl^db) = ^(B)^), = y|cosh(|(z/2 - £))r(i£ + i) 

Essentially S is the S-matrix, and TZ(E) is the reflection coefficient (phase shift). 



(3.6) 



The closed string momentum modes correspond to operators in MQM of the form[43 



V ±k/R (x) = e ±ikx / R TrX±(x) k / R (3.7) 

where x is the Euclidean time variable, R is the radius of the thermal circle, and the 
allowed momenta are ±k/R for integer k. A general perturbation of momentum modes is 
described by a potential 

V±(x±) = Rj2t±kX k J R (3.8) 

fc>i 

In the c = 1 MQM deformed by (|3.8|), the energy eigenstates are given by "dressed" wave 
functions 

*±(s±) = e^^'^tyf (z±) = W±Vf (z±) (3.9) 
where the phases ip± are of the form 

<P±{x±;E) = V±(x±) + \(t>{E) -Rj2lv ±k (E)xl k/R (3.10) 



k 
k>i 
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In above <p(E) is a constant phase shift. The terms involving v±k vanish as x± — ■> oo. 
Semiclassically the perturbation (|3.8|) corresponds to deformed fermi sea profile 



x -k/R (3.11) 



x + X- = n + x±d(p±(x±; /j,) 

= y- + ^2 kt ±k x ± R + ^2 v±k ( E ) 



k>l k>\ 



The compatibility of the two equations with + and — signs in the subscripts puts con- 
straints on v±k in terms of t±f e , s. This comes from the requirement that are the wave 
functions of the same state in x± -representations. They are related by = or 

equivalently 

W_ = W+ti (3.12) 

Note that this constraint is essentially determined from the reflection coefficient 71(E). 

To see the integrable structure of Toda lattice hierarchy, we shall recast the above 
in terms of operators on the -E-space. For example, it follows from ( |3.5|) that x± are 
represented as shift operators a)^ 1 = e^ ldis . One can define a Lax pair 



V fc >* 



(3.13) 



k>l k>l 



which satisfy the commutation relation 

[L ± ,M±} = ±iL± (3.14) 
Recall that the dressing operators W± in terms of E and Cj are of the form 

W± =e^ 2 ^l + E w ±^ Tfe/ii j e T ^^>^ ±fe " ±fe/ii (3.15) 

The integrable flow equation of the Lax operators is 

dt n L± = [H n ,L±] (3.16) 



where 



H n = (d tn W+)W+ l = (d tn W.)WZ x 
11 



(3.17) 



By the virtue of ( |3.12|) , the generator of the Toda flow H n is the same for both L + and 
L-. By a standard argument (for example, see p3|), it follows from the structure of ( |3.15| ) 
and ( |3.13| ) that H n are of the upper or lower triangular form 



"<\ _i_I<T'"'" 1() . 

n > 0) (3.18) 



H n = (L^ H ) > + -(L + ,„ 



H- n = {L n J R ) < + 1 -{L n J%. 

It also follows from (|3.16|) that the zero- curvature conditions on H n are automatically 
satisfied. 

The "almost lower triangular" structure of lax operators L± (|3.13|) , together with 
( |3.160 and the important relation (|3.18| ), define a Toda lattice hierarchy. The two Lax 
operators L + [t n ] and L-[t n ] are in addition constrained by ( |3.12Q . Using the functional 
relation R(E — i) = (-E + i/2)R(E), the constraints on Lax operators can be expressed 
in terms of the so called string equations, in this case given by 

M+ = M_ = -(L + L_ + L_L+), 

+ 2 l + +h (3.19) 

[L+,L_} = -i. 

These conditions define a constrained Toda hierarchy. 

The parameters 0, v±k appearing in W± are determined in terms of t±k via the Toda 
flow ( ETT5D and the initial condition L±\ tn =o = e T ^°/ 2 c2; ±1 e ± ^ / 2 where e ^ o(E) = 71(E). 
They are related to the Toda r-function r(i?;i n )i by 

Vn = S jEl, m = i^ E + ^ R \ (3.20) 
dt n ' vy ' r(E-i/2R) K ! 

where we suppressed the dependence on t n 's. The r-function satisfies Hirota's bilinear 
equations [[16], 44], an infinite set of differential equations in £ n 's that completely determine 
r(E; t n ) provided the initial condition r\t n =o(E) for all E. The initial condition for r(E; t n ) 
is essentially equivalent to the constraints imposed by the string equations. 

To see how the r-function is related to the free energy, let us compute the density 
of states in terms of <j>(E). We shall introduce a cutoff at x = a/2A, which is a wall 



The r-function appearing here is denoted r'[t] in [44]. In section 3.4 we will introduce 
another r-function r[t] defined using vertex operators, following [45,21]. They are related by 
r'[t] = ex P(Xln>o n ^™^-™) ' r M- The grand canonical partition function is given by r'[t]. Also the 
standard form of the r-function, n[t], is related to T(E;t n ) by n[t] = r(/i + il/R;t n ). 
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that reflects all momenta. In light cone variables this imposes a boundary condition at 

£ + = X- = a/A, 

9*(VX) = 9*(VX) (3.21) 
as defined in ( |3.9|) have asymptotic behavior 

*|(a/A) ~ e Ti ^/ 2 (VI) ±iE x (£ independent piece) (3.22) 



It follows form (|3.21|) that the allowed energies E n satisfy i(p(E n ) — iE n In A = lit in. The 



density of energy eigenstates p(E) is then 



p(E) = H± _ ±.*m (3.23) 



The free energy is given by 

/oo 
dEp(E) ln(l + e- 27ri?(E+ ^) 
-oo 

<j>(E) 



-oo 

"OO 



R / dE- 



1 _|_ e 2irR(E+/j,) 

^(t){-a + {n+\)i/R) 



(3.24) 



n>0 



= lnr(-/i;£ n ) 

where in the third line we closed the contour in the upper half plane which picks up the 
poles at E = —u + (n + ^)i/R, and in the last line we used the relation (|3.20| ). In our 
convention a is the negative chemical potential, which is positive when the fermi level is 
below the top of the potential (E = 0). We see that the perturbed grand partition function 
is Zp[t n ] = r(-u;t n ). 

3.2. "Theory II" 



Now let us turn to type OB MQM, also known as "theory of type II" ||47|| . For the 



unperturbed Hamiltonian there are two sets of eigenfunctions that classically correspond 
to fermions in the left and right sector of the fermi sea 

i 



±iE 

1>+>(x±) = —== X± =, (x± > 0) 



(3.25) 



2tt a/1 + e 27rE 
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In the other half of the real x± axis the wave functions are defined by analytic continuation, 
explicitly 

il>l > {-x ± ) = ±ie* E ij>l > {x ± ), 

(3.26) 

^±,<(z±) = ±te* E i;l < (-x ± ), x± > 0. 
The wave functions in x+ and X- representations, if)+ and i/j-, are not merely related by a 
reflection because of the quantum tunnelling. The operator S relating to ijj- is defined 
as ( |3.3| ) but with a different kernel 

K(x-,x+) = -—e ix + x - (3.27) 



2,7T 



S acts on the energy eigenstates as [23 



where 

= ^Le-* (E - l / 2) r(z£ + i) (3.29) 
V 2tv 

Perturbatively the left and right side of the fermi sea decouple. This is reflected in the 
exponential suppression factor e KE (E < 0) when we analytically continue x± — > —x±. S 
is diagonalized by symmetric and antisymmetric eigenfunctions 

, f „ = ^±^, ^, =fc^ (3.30, 

Then we have 

SiP e >s (x_) = R 8 (E)ip E >s (x_), R S (E) = y|cosh[^(z/2 - E)\Y{iE + \ 

$V+,a(*-) = i?a W- |0 (flJ-), = ^mnh[^(i/2 - E)]T(iE + §). 

Note that x± act on ip E > as cD ±:L = e^ ldE and on < as — u^ 1 = —e^ ldE , or equivalently 

a±^f,. = V±?, *±^f,«=^±? (3-32) 

Now let us turn on perturbations by closed string momentum modes. The dressed 
wave functions are of the form 

*±>±) = e^^ E ^l s (x ± ) = W&l.[x ± ), 

14 



2)i 



(3.31) 



where 

k>l k>l 

~ (3.34) 

<PW^ X \E) =Rj2t a ± k " ±k/R + ^ a {E) -Rj2 % v a ±k (E)u Tk/R 
k>i k>i 

t± k ,t± k parameterize the symmetric and antisymmetric part of the asymptotic pertur- 
bations of the fermi sea, respectively. W± again have the structure of ( |3.15| ). The Lax 
pairs (L±, M±) are defined the same way as ( |3.13|) , and satisfy the relations ( |3.14| ), (|3.16|) , 
( P-17|) . Now we have two sets of generators of the flow, and (associated to and 
t®). Since L± and Cj are block diagonal in the (ip s ,ijj a ) basis, 

(K \ 

H n = 3.35 

V o K) 

still satisfy the upper or lower triangular relations ( |3.18| ). They define two independent 
Toda flows, associated to the symmetric and antisymmetric perturbations of the fermi sea. 

However, cD ±1 is no longer the same as x± since the latter interchanges ip s with ip a . 
Consequently the functional constraint on the reflection coefficients is modified to 

R s (E-i) = (-E + i/2)R a {E), 

(3.36) 

R a (E-i) = (-E + i/2)R s (E). 

So the strings equations (|3.19|) no longer hold in the case of "theory II" . They are modified 
to 

M s,a = M «,a 

L s + L s _ = tanh[^(z/2 - M s ± )] (— MJ + i/2), 

L S _L S + = coth[|(i/2 - M±)] (—M± - i/2), (3 . 37) 
L\L a _ = coth[-(z/2 - M£)] (-M| + i/2), 
L a _L\ = tanh[|(z/2 - M£)] (-M| - i/2). 

Perturbatively, i.e. in the limit e nE <C 1, these equations reduce to (|3.19|) . They define 
two independent constrained Toda hierarchies, with the Lax operators acting on and 
ijj a respectively. The perturbed grand canonical partition function is the product of the 
two r-functions, 

Z t ,[t s ,t a } = T s (fx;t s )r a (fJ,;t a ) (3.38) 

As remarked in [E5|, this doesn't mean that the NS-NS and R-R closed string modes 
decouple from each other. The symmetric and antisymmetric perturbations of the fermi sea 
are mixtures of NS-NS and R-R fields in spacetime, due to the nonlinearity of bosonization. 
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3.3. Type OA theory 

Type OA MQM can be represented by non-relativistic free fermions moving in a two 
dimensional upside-down harmonic oscillator potential. The Hamiltonian is 

H=l(Pl+P 2 y )-^(x 2 + y 2 ) (3.39) 

The theory has different independent sectors labelled by net DO-brane charge q, which is 
the same as the angular momentum J = xp y — yp x || . We shall mostly focus on the case 
where there is no net DO-brane charge, namely the J — sector. 
It is again convenient to define light cone variables 

J—r X ± p x —F=rV ± Pv 

X±= V2 ' V±= V2 (3 ' 40) 



and 



We have commutators 



z±=x±+iy±, z± = x±-iy± (3.41) 



[x+,x-] = [y+,y~] = i/V2a', 

[z+,z-] = [!+,§_] =0, (3.42) 
[z+, Z-] = [f+, S-] = 2z/v / 2c7. 
In light cone variables the Hamiltonian is written as 

H = -(x + x- +y+y- - 



V2a> 

+ ~ -7=) (3.43) 



2^ ^ 
i , d d 1S 

= 7l7 (z+ a^ + " + si7 + l) 

where the last line is written in (z+, ^representation. In addition, we have commutation 
relations 

[J, z±] = z±, [J, l±] = -l±, (3.44) 

or in (z+, z + ) representation 

j =*+£- + -* + £- + (3 - 45) 

The wave function of a state can be expressed either in (z + , z + ) representation or in (z_, z—) 
representation, denoted by ifj + (z + , z + ) and (£_,£_) respectively. Since we restrict our 
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wave functions to have zero angular momentum, it must be of the form ip±(z±,z±) = 
ifi±(z±z±). The energy eigenstates are given by (setting a' = 2) 

, 77 HE — i zniE— i , N 

</>f~2± 2 z ± 2 (3.46) 
The wave functions in (z+, z + ) and (z-,Z-) representations are related by 

if)-(z-,z-) = (Sip + )(z-,z-) 

= J dz + dz + K(z- 1 z + )K(z- 1 z + )ip + (z +1 z + ) 

where K(z-,z + ) = -^=e lz ~ z+ . Acting on energy eigenstates, we have 



(3.47) 



■fy? = K(E)^ K(E) = r, {iE *jL <•>• l8 ) 

1 \—irLi + 5 J 



This is the same as the phase shift found in H. 

Now we can study perturbations of the fermi sea. Let us consider dressed wave 
functions 

where the phases ip± have Laurent expansion 

<p±(z±z±;E) = \^{E) + RY J t±k{z ± z ± ) k ' R ~ RY J \v±k{z ± z ± )- k ' R (3.50) 

k>l k>l 

t±k parameterize the asymptotic perturbation by momentum modes of NS-NS scalars, 
corresponding to the operator 

V p = k/R = Tr(Z ± Z ± )W R (3.51) 

in type OA MQM, where the sign of the subscripts depends on the sign of k. (z±z±) can be 
represented as shift operators a) ±1 , where u acts on energy eigenstates as a) ±1, i/'± = ^±^ 1 - 
We have cD = e~ ldE and the commutation relations 

[u) ± ,-E]=±iu ± (3.52) 

As before we can define a Lax pair 



L± = W±u ± W± 1 , M± = -W± EWg 1 
17 



(3.53) 



The dressing operators W± satisfy the constraint 



W- = W+ ■ 71(E) (3.54) 

This again defines the structure of constrained Toda lattice hierarchy. But the string 
equation Q3.19|) is modified to the following 

M + =M_, [L+,L_] = 2zM ± , {L + ,L_} = 2M| - -. (3.55) 

The density of states, and hence the free energy, is related to the phase <p(E) in the 
standard way. To compute the density of states, we shall introduce a cutoff at x 2 + y 2 = A. 
The cutoff wall reflects all the momenta, so we have xp x + yp y = as well. Further we 
demand the vanishing of angular momentum xp y — yp x = 0. The combination of these is 
equivalent to z + z + = z-Z- = A. We can impose a boundary condition at the wall 

#+(A) = tf^(A) (3.56) 

It follows that the density of states depends on E the same way as in ( |3.23| ). 

It is not hard to generalize the above construction to sectors of nonzero net DO-brane 
charge q. These backgrounds are identified as extremal black holes in type OA string 
theory [|12"|| . In this case, the energy eigenstates (carrying angular momentum q) are 



E ±iE+q/2-\ _±%E-q/2-\ 

y)± ~ z± z± (3.57) 
The reflection coefficients are computed from ( |3.47| ) to be 

= [ \ (3.58) 

The Lax operators are defined as before, but the string equations become 

a 2 - 1 

M + = M_, [L + ,L_] = 2zM ± , {L + ,L_} = 2M| + . (3.59) 

Interestingly, {L + ,L_,M = M±} form a representation of sZ(2,R), with isospin I = 
(M-l)/2- 
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3. 4- Type OB IMM perturbed by closed string momentum modes 

Following ]2y,^TJ, we shall integrate out the tachyons t( n \ t^ in the type OB IMM 
( p.lO| ) by analytic continuation, and obtain a matrix integral in terms of the eigenvalues 
of X and X, 

_ _ OO i 

/>i . nn . sy . ry . j 



n^n-n 2 (^)sm 2 ( ^)n n 



(2nnR + Xi — Xj) 2 — a 2 

i<j i, j n=—oo v J/ 

= [ TldxdxTlsin 2 ( Xl ~ Xj )sin 2 ( Xl ~ Xj )T\ - 

J 2R 2R sin [(xi - Xj)/2R] - sin 2 (a/2 J R) 

(3.60) 

Roughly speaking we have a system of iV eigenvalues Xi and Xj , they repel eigenvalues of the 
same type through the F-P determinant and interact through some effective potential. The 
analytic continuation made above is quite naive. As a trade off, ( |3.60| ) is not unambiguously 
defined, and we have to give a correct contour prescription. This will be considered in the 
next section, for now we still formally work with ( [3.60|) . 

Let us write the partition function in terms of the integral over eigenvalues of U, U, 

/y _ f^^j I 'y ■ (3 ^ ^ j I 



N 

1 y-r I dz k I dz k y-r 



I %i Zj | | Zi Zj | 



%3 

N 



Ziq 1 / 2 — Zjq 1 / 2 \ ■ \ziq 1 / 2 — Zjq 1 / 2 ] 



(3.61) 



v > k=i" " i^j 

U 



4-4 (ziq 1 / 2 — zjq l / 2 )(ziq 1 / 2 — Zjq 1 / 2 ) 

where q = e %a / R = e 7r «v / 2a 7 /-R Consider the perturbation by momentum modes (winding 
modes in the T-dual type OA theory) 

AnTrtT* + ~X n TrU n (3.62) 

n6Z 
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The generating functional is then 

ma, a] = * n / S / %e«»»™ n (* - *x* - %> 



x n * 

11 (Z^V2 _ ^-l/2)(^-l/2 _ ^ 1/2) 
1 rr /" g^i. f dz 



(3.63) 



n 



(iV!) 2 / 27T2 7 2ttz 



x det 71 I — -t^ — T-77T ) det 



where 



u( 



i(z) = J2 X nZ n , u{z) = Y,~K~z n , (3.64) 

n n 

and we have used the Cauchy identity 

= det , 3 ( 1 ) (3 . 65 ) 
It is most convenient to consider the grand canonical partition function 

oo 

Z„[A, A] = J2 e wV ^^ N Z N [X, A] (3.66) 

JV=0 

where /U is the chemical potential of type OB theory, related to the one of the T-dual type 
OA theory by /j,' = fj,R/ \j2ot' . For reasons that will become clear shortly, let us define 

A n = 2t n + (q n + q- n )t n , ~X n = -2t n - (q n + q~ n )t n (3.67) 

and a "r-function" 

oo 

Ti\t t] = e~ n[2t„t_„+2t„t_„ + (g"+ 9 -")(t„t_„+t„t_„)] ^21 ^V^oF^N [t t] 

n=o ' (3.68) 

_ -Y\ n[2t„t_„+2t„t_„ + (g"+g-")(t„t_„+t„t_„)] 7 r, 7] 

One might hope that the grand canonical partition function [t, t] is the r-function of 
some integrable hierarchy. This is not quite true. The perturbations that generate the 
integrable flows are related to £ n ,£ n 's through some bosonization maps, as we will show 
below. 
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It is useful to rewrite the partition function in vertex operator formalism, analogous 
to the case of c = 1 string [21]. One can introduce two independent 2D chiral bosons 
<Pi,2{ z )-i with mode expansion 

o-(l,2) 

<Pi,2(z) = qi,2 + Pi,2 lnz + J2 ~^^z~ n (3.69) 

The vacuum \l) is defined by 

^ 1} ' (2) I0 = (n>0), p 1 \l)=p 2 \l) = l\l). (3.70) 

We could have considered more general vacuum state l 2 ), but that would be unnecessary 
for our purpose. We further define 

<Kz) = <p 1 tf' 2 z)-<p 2 (q- 1 / 2 z), 

3 - 71 

<t>(z)=^(q- 1 / 2 z)-^ 2 (q 1 / 2 z), 
and H n , H n the corresponding creation and annihilation operators in the mode expansion 

Of 0, 0, 

H n = q- n ' 2 H^ - q n ' 2 H^ 
H n = q n ' 2 H^ - q'^H^ 



(3.72) 



Note that 4> an d 4> are n °t independent fields. We have commutation relations 

[H n , H m ] = 2n5 n+m = [H n , H m ], [H n , H m ] = (q n + q- n )n5 n+m . (3.73) 
Using the operator 

the "r-function" ( |3.68|) can be written as 

n[t,t] = (/|e"^»>o t ™ H " + *™ H ™g / e E„>o t -™ H -™+ t "-" H -"|Z) (3.75) 

To see this, observe that the contractions among the operators and in (|3.74| ) give the 
integrand in ( |3.61[ ), and commuting g' through e~ ^™>° i ™ Hn+i ~ ni ^™ and e^™<o inH " +i ~"^ n 



give the momentum mode perturbation e u+u appearing in ( |3.63|) . The symmetry (p — > <p+a, 
4> — > 4> + a is respected by the vacuum state [/). Therefore, we can replace g' by 

g = exp (q-W/ 2 j ^- : e *(*) + e~*W :^ (3.76) 
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in the partition function (|3.75|) . The partition function in sectors of nonzero net D- 
instanton number can be obtained from the general vacuum state | Zi , Z2) - The latter 
are closely related to the solitonic sectors of the type OB bosonic Hilbert space discussed 
in |25[| . Of course, the r- function ( |3.75| ) does not simply factorize into two components 
involving only t n and t n respectively, because <p an d <fi have nontrivial OPE. 

The exponent in g can be interpreted as a Hamiltonian. It is clearer to fermionize 



e Vl(2) ~ ^i(z), e"^ l(z) 



(3.77) 



The vacuum |/) satisfies 

ipx,r\l) = 4,-r\l) = M l ) = 4,- r \l) =0, r e Z + §, r > I. (3.78) 
In terms of i/ji,2, we can write 

'-^{q-^z^iq^z) + MQ~ 1/2 z)MQ 1/2 z)] } (3-79) 
We can interpret (|3.75|) as a partition function of the fermions V 7 2- The integral in the 



exp < q l ^ R / 2 



dz 
2^ 



exponent of g can be written in first-quantized form as 

_ g| In q(zp z +p z z) p 

(3.80) 

= e -Tk( z Pz+P* z ^ P 12 

where P12 : ipi — > ip2i ^2 — ^lj and we have used g = e %a l R . This "Hamiltonian" can be 



compared to the c = 1 string case|2T||, where one can write the partition function in terms 
of vertex operators involving a single chiral boson ip(z) or its fermionization i()(z), and the 
corresponding "Hamiltonian" is H = e~^^ zpz+PzZ K 

Since ip± = 775 C0i ± ^2) diagonalize P12, it is clear that Ti[t n = t n = 0] factorizes as 
the product of partition functions involving ^_|_ and ijj- separately. In the type OB fermi 
sea picture, should correspond to symmetric and antisymmetric perturbations of 

the fermi sea. Indeed, as we will show in the next section, the partition function involving 
^4. (ip-) agrees with the partition functions in type OB MQM involving only symmetric 
(antisymmetric) perturbations. 

In the Hamiltonian interpretation, the "initial" state in the expression ( |3.75|) for the 
r-function is 
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where 

4 1} = Q~ n/2 tn + q n/2 t n , ti 2) = -q~ n/2 t n - q n ' 2 t n (3.82) 

As we have seen, ifii ± ifa correspond to symmetric and antisymmetric perturbations of 
the fermi sea. It is natural to expect that -01 and %j)2 correspond to perturbations of the left 
and right sector of the fermi sea (decoupled perturbatively) . Consequently the bosonized 
modes Hn\ are linear combinations of momentum modes of NS-NS and R-R scalars 
in spacetime. From ( |3.67| ) and ( |3.82|) , the perturbation parameters <4 , *n are related to 
X n , X n in IMM by 

An = q n/2 tP - q-/HW 
An = -q- n/2 t^ + q */HW 
This also leads us to the identification 



(3.83) 



a L ,n ~ q n / 2 Tre inX/R - q- n / 2 Tre inji/R , 
a Rn ~ - q - n / 2 Tve inX/R + q n / 2 Tre inX / R , 



(3.84) 



or 

q 1 a,L, n + q 1 clr,' 



(3.85) 



q n _ q -n 

r^jnX/R _ q~ n/2 aL,n+q n/2 aR,r 

~ q n _ q -n 

where aL, n ± a>R,n are the momentum modes of the NS-NS and R-R scalar with momenta 
p = n/R. Note that the Lorentzian energy is E = in/R. To analytically continue to 
Lorentzian signature, we have q n = e 2nE . Perturbatively positive powers of q can be 
neglected (E < below the top of the potential). In this limit we have approximately 
rj^^nX/R ^ —e wE a,R in , r £re inX ^ R ~ —e 7TE a,L,n- As will be shown in section 5, the relations 
( |3.84[ ) precisely reproduce the ZZ boundary state in super-Liouville theory. 

We should remind the reader that i are not the bosonization of ± iip%, 

and that tn ± ttn^ are not the same as t n ,t n defined in section 3.2. The operators 
r £ re inX / R ^ r^ e inX/ R - m iy[y[ correspond to linear combinations of NS-NS and R-R modes, 
whereas in MQM the fluctuations of the fermi sea are related to closed string fields in 
spacetime by bosonization. Due to this bosonization relation between t n , t n and t n ,t n , the 
partition function of IMM does not factorize in a manifest way. However if we replace the 
"initial state" (and similarly the "final state") in (|3.75|) by 



"10 (3.86) 
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where are the modes of the bosonization of ip±, then the partition function will factorize 
into two r-f unctions of constrained To da hierarchies that depend only on t s or t a , just as 
in type OB MQM. To show that the perturbed grand partition functions in type OB IMM 
and MQM are actually the same, it remains to show that the constraints (string equations) 
of the IMM agree with those of MQM. 

3. 5. String equations for IMM 

In this subsection we derive the string equations for type OA and OB IMM, by com- 
puting the commutator [L+, L_] at t n = 0. Using the general form of Lax operators ( |3.13| ), 
we know that in the absence of perturbation, 

L± = e ^o/2 c2 ,±i e ±^o/2 ) 

(3.87) 

[L+, L-] = e i( t>o{E+i)-i<t>o{E) _ e i4> {E)-i<t> {E-i) ^ K ' 



Using the relation ( |3.24| ), we have 



(3.88) 



where T(fi) = \n.Z(fi) is the unperturbed free energy. To show that the IMM string 
equations are the same as the ones for MQM found in previous sections, it suffices to 
show that 0(/i) as given by (|3.88|) (with T evaluated from IMM) does lead to the correct 
reflection coefficient 11(E) = e ^ (E) . 



Let us first consider c = 1 or type OA IMM, which is already done in []20j . The grand 
partition function takes the form of a Freholm determinant 

Z([i) = det(l + q iliR k) (3.89) 

where q = e 27ri / R , \i is the chemical potential of the theory on radius R (the chemical 
potential of the T-dual theory is // = fxR). K is defined by 

(Kf)(z) =-<f ^ fiZ \ 2 , (3.90) 

The contour prescription of [^DJ (see also |41|]) is to add a small imaginery part to R so 
that |g| < 1. A basis that diagonalizes K consists of the monomials z n , with 

, q n+ iz n , n>0 
Kz n = { ~ (3.91) 

I, n < 
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The free energy is given by 



Jfy) = ^ln(l + ^ i? + n +2) 



n>0 

oo 



' 771 q m /^ — g TO/2 



(3.92) 



/ c t 4sinh(t/2)sinh(t/2.R) 
where C is the contour that picks up the poles of sinh(t/2) at t = l-nni (n > 0). For 
example, we can choose C to run from — oo to along the real axis, and from to zoo on 
the right of the imaginery axis, as shown in Fig. 2. 



Fig. 2: The integration contour C in the complex i-plane that appears in 
where R is given a small negative imaginery part. 

( p.92|) is however not quite right. It is in fact the free energy of a system with ordinary 
("right sign") harmonic oscillator potential. Nevertheless, the integrand in the last line of 
( fj.92| ) has the same structure as in the free energy computed from the true spectrum of 
the upside-down harmonic oscillator potential. This suggests that we modify the contour 
prescription, so that 

T = Re — —-^ = , ; (3.93) 

J t 2sinh(t/2)sinh(t/2i2) V ! 



It then follows from (|3.88|) that 

1 dt sin(/i£) 



o t sinh(£/2) 
£g/x+j) 

r(-v+l) 



(3.94) 

-i In 
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which agrees with the type OA MQM reflection coefficient (|3.48|) . 

The change of contour in above calculation, which effectively takes us from the right 
sign harmonic oscillator potential to the wrong sign potential, seems rather ad hoc. It 
will be nice to have a clearer prescription. We will however use the same prescription to 
compute the free energy for type OB IMM, and find a nontrivial agreement with the results 
obtained from the MQM. 

In the case of type OB IMM, we have seen from the previous subsection that the 
unperturbed partition function is the product of two Fredholm determinants, 

Z{ji) = det(l + iq if * R / 2 K) det(l - iq lflR/2 K) (3.95) 

where again, \x' = fiR is the chemical potential of the T-dual type OA MQM. The "sym- 
metric part" and "antisymmetric part" of the free energy are 

F a (y) = ln(l + iq ^ R / 2 + n +\) 

n>0 

dt e iM*/2 e t/4 



! c t 4sinh(t/2)sinh(t/2i?)' 

n>0 

dt e *M*/2 e 3*/4 



c t 4sinh(i/2)sinh(t/2i?)' 

where C is the contour that picks up the poles of sinh(t/2) at 2nni (n > 0). Again, with 
proper modification of the contour prescription, and using ( J3.88[ ), we obtain the phase 
shifts 

dt sin(^t/2)e t / 4 
o T sinh(t/2) 



(f>s (m) = 

= — i In 

<f>a(V>) = ~ 

= —i In 



2 4 (3 97) 

° dt sin(^/2)e 3f / 4 

t sinh(t/2) 

r(f + |) 



n-f + !) 

Up to terms that contribute a constant to the density of states, these precisely reproduce 
the reflection coefficients ( |3.31|) . 
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Note that one can assign different chemicai potentials fj,± to the symmetric and anti- 



symmetric perturbations. This corresponds to turning on a background RR flux [[16]]. The 
T-dual type OA flux background should be analogously defined by giving different chemical 
potentials to two independent winding sectors. In this case we need to work with the grand 
canonical ensemble which sums up sectors of arbitrary DO-brane and anti-DO-brane num- 
bers. We don't expect the corresponding type OA background to have a fixed net DO-brane 
charge q. The thermal fluctuation in q should be of order (9(e _7r ^ M ) instead of 0(e _27r ^ M ), 
since the "effective" chemical potential for the two winding sectors in type OA theory on 
radius 1/R is /j,R/2 (// being the chemical potential in the T-dual type OB theory), as in 
( p.95|) . This is responsible for the mismatch at nonperturbative level between type OA 
partition function with (j,± = [i ± Q and type OB partition function with RR flux q = iQ 



as found in [16 



With proper contour prescription for the type OA and OB IMM, we have produced the 
correct reflection coefficients, hence the string equations of the constrained Toda lattice 
hierarchies. This proves the equivalence between type OA (OB) MQM and IMM perturbed 
by purely momentum modes or winding modes. The contour prescriptions introduced 
above should be regarded as part of the definition of the double scaling limit for the IMM. 



4. Connecting IMM to MQM via Tachyon Condensation 

4-1. The array of D-instantons 

In this section we shall attempt to connect the IMM to MQM in a more direct and 
intuitive way via open string tachyon condensation. To be definite let us work with type 
OB string theory compactified on a thermal circle of radius R = ma, where m is an integer, 
a = iry/2a' is the critical distance. As well known |26|, turning on an open string tachyon 
profile 

T(X) = Xcos(nX/a) (4.1) 

on the Euclidean world volume of an unstable DO-brane is described by an exactly marginal 
deformation in the boundary CFT. For sufficiently large A one ends up with an array of 
alternating D- and anti-D-instantons separated at distance a. From the point of view of 
matrix models, this suggests that the MQM perturbed strongly by the tachyon profile 
( |4.1| ) should be identical to the IMM expanded near the configuration of the array of D- 
instantons. Although it is not clear to us how to describe the tachyon profile (|4.1| ) in the 
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MQM in a precise way, it is very plausible that the effect of ( f4.1|) is to put the MQM on a 
discrete Euclidean time lattice of spacing a.i 

Given this picture, we want to understand how the operators on both sides are identi- 
fied. When the open string tachyon is condensed, the DO-brane effectively vanishes, so the 
open string excitations are localized near the "sites" of the time lattice where T(X) ~ 0. 
Let us consider in the MQM a small tachyon lump ~ near x = (k+\)a. This will 
shift the zero locus of T(X) = \cos(%X/a) + to x k ~ (k + |)o + c(-l) fc $ fc , where 

c = a/irX. So effectively the position of the kth D-instanton (or anti-instanton) is shifted 
by c(— l) k <&k- This suggests that the positions of the D-instantons in the array in IMM 
should be mapped to the open string tachyons on the corresponding sites in the discretized 
MQM. Let us denote by Xk the fluctuation of the position of the kth cluster of instantons 
from the array configuration. Then we are tempted to identify 

X k ~ (-l) fc $ fc (4.2) 

On the IMM side we shall integrate out the (complex) tachyons which are open strings 
stretched between D- and anti-D-instantons, and get an effective theory that describes only 
the collective coordinates X k of the instantons. The D-instanton and anti-D-instanton want 
to move toward each other so that the tachyon can condense. Roughly one can think of 
this instability as an unstable "effective potential''^ V{X) (of periodicity a) felt by the 
eigenvalues of X^s. Similar to the picture of c = 1 MQM, here the eigenvalues repel each 
other and fill up the "valleys" of the potential. In the large iV limit, the eigenvalues are 
distributed in 2m cuts along the thermal circle. When there are a sufficiently large number 
of eigenvalues, the cuts will connect to adjacent ones and we expect a phase transition. The 
transition point is where the "energy" of the eigenvalues reaches the top of the potential. 
From ( |4.2| ) we expect that this is the critical point that defines the double scaling limit! 
In other words, the array of D-instantons at critical distance in IMM plays the same role 
as the unstable DO-brane corresponding to the top of the tachyon potential in MQM. 



3 It is well known [|TJ that the discretized MQM on a time lattice of spacing e < a is exactly 
equivalent to the continuum theory up to a redefinition of parameters. In the critical limit e = a 
the discretized MQM becomes singular. We will regularize the theory by taking e slightly less 
than a. 

4 The picture we are describing here is rather heuristic. As we will see in the next subsection, 
the interaction between the eigenvalues are not simply represented by an ordinary potential plus 
the repulsion through the Vandermonde determinant. The phase transition occurring here as 
the eigenvalues merge the top of the "potential", is similar but different from the Gross- Witten 



transition [48], the latter being well known to describe 2D gravity coupled to c = matter|2 



28 



4-2. Effective matrix integral of the array 



We will expand type OB IMM near the configuration of an array of alternating D- and 
anti-D-instanton clusters at separation e. It is natural to work with e > a so that all the 
tachyons stretched between D- and anti-D-instantons are massive. This can be thought of 
as a kind of regularization. Eventually we will be interested in the limit e — > a. We will 
assume the radius R = 2ma/2n, so there are 2 m clusters distributed along the thermal 
circle, with N instantons in each cluster. We shall write the diagonalized X and X as 



X = 



\ 



X 2 + 2e 



\ 



X = 



X m + (2m-2)eJ 



\ 



X 2 +3e 



X m + (2m - l)e / 

(4.3) 

where X^ = diag{x2fc-i, a }^ = i and X^ = dia,g{x2k,a}a=i are the collective coordinates of 
each cluster of N instantons. 

We will take e to be very close to a, and ignore all the massive modes with masses of 
order ~ a, so that only the "tachyons" as open strings stretched between adjacent clusters 
are retained. Furthermore, we shall take the double scaling limit before taking the limit 
e — > a, so we can assume the fluctuations x^ a are small, since only the small fluctuations 
near the top of the potential are responsible for the universal behavior of the theory near 
the critical point. With these approximations, the (unperturbed) partition function can 
be written as 

9 , ii _^ 



J N 



where Zj a = e l( -- Xj ' a+ ^ e ^ R : q = e ia l R . We can rewrite it as 



q 1 / 2 z j+1 ^ - q 1 / 2 Zj jl 3 



(4.4) 



\\ dx jiDl det Q/3 exp 

j,ot j=l 

/2m 
Y[dx ja Y[ det Q/3 exp 



lnsin(— — ; ) 



2R 



1 



(4.5) 



Let us note that there is a zero mode corresponding to the overall shift of X and X, whose 
origin can be traced back to the decoupled diagonal £7(1) in the U(N) x U(N) type OA 
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MQM. We are free to fix this redundant gauge symmetry. Then we expect a tachyonic 
instability for each individual Xi and Xj as they want to move towards each other. 
Using the Itzykson-Zuber formula, we can rewrite ( |4.5| ) in the limit e — > a asi 

2m ^ 

!l ^£.T*(x, +1 -o,x j n+ ) 2 



3=1 



/tilth 
\ \ dx J dn J 5(x J+1 - n^nt) 

3 = 1 



(4.6) 



In fact, by integrating out all but one of the X,-'s, ( f4.6| ) reduces to the three matrix model 

of h 

1 rfX + dX_rfOe lTr ^+ x --^+ nx - nt ) (4.7) 

with q = e 27 ™ fl /v / 2o 7 Q ur p rescr jption is to compute correlators with imaginery R, and 
then analytically continue to real values of R(— m\j2a' in above). 

4.3. The discretized MQM 

Now we put the MQM on a discretized time lattice of spacing e, where e < a. In 
order for the discretized MQM to be exactly equivalent to the continuum MQM, we need 
to have discretized propagator 



dU iii+1 exp -p'eTr 



U^ +1 $ l+1 U t / +1 - $ 



\2 1 , J2(jr.2 



(4.f 



where dU^i+i is the Haar measure over U(N), and the parameters p',u>' are related to 
those of the continuum model by 



we = 2sm— , p = p~— — 

2a sm(7re/a) 



(4.9) 



We can diagonalize $^ into its eigenvalues Xi^ a at each step, so the propagator simplifies 
to 



dll exp 



0' 



a a, (3 a. 

(4.10)' 



5 To compare with the three matrix model, we are really taking e to differ from a by a small 
imaginery number. 
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Using the Itzykson-Zuber formula 



/ 

Ju 



N-l 



'U(N) 

The propagator becomes 



A(x)A(y) 



(4.11) 



k=i 



exp -^(l-|o;'V)^ +1>a + At; 



A(A l+1 )A(A,) 



(4.12) 



The partition function is given by 

2m 



J n rfA ^li de Wexp{-^ [(A i+lja -A^) 2 -|o;' 2 6 2 (A? +l5a + A^)]| (4.13) 



Note that we haven't made any approximation in above manipulations. Under the identi- 
fication (gT9j), (g7L3|) is exactiy the same as the path integral of the continuum MQM. 
In the limit e ~ a, ( 4.1 3| ) is approximately 



2m 



J\ d Ka ]^[ det Q/3 exp 



,3 



a — t 



-(A m , Q + A lj/3 ) 2 + ^_(a - e) 2 (A 2 +ljQ + A 



(4.14) 

Under the identification (4.2) , (|4.14| ) approaches the path integral of IMM expanded near 
the array (33), up to rescalings of e and A^'s which can be absorbed into a shift in the 
chemical potential. This confirms the connection between the discretized MQM and the 
IMM expanded around the D-instanton array in the critical limit. 
Similar to the manipulation of ( [4. 6| ) , we can rewrite ( [4. 14j ) as 



JJdXjdYjdQje 

3=1 



(4.15) 



which again reduces to the three matrix model. The contour prescription of giving R a 



small imaginery part so that |g| 
imaginery number. 



< 1, corresponds to taking e — a to be a small 



5. On SD-branes and ZZ Boundary States 



It was pointed out in []28| that sD-branes (A = 1/2 s-brane) can be described as 
an array of D-instantons along Euclidean time separated at the critical distance. In the 
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previous section we have argued that sD-branes play the same role in the IMM as D0- 
branes in MQM. In this section we study the closed string fields sourced by sD-branes 
from IMM. We find agreement with calculations from ZZ boundary states. This is very 
much in the same spirit as the calculation of D-brane decay into closed strings from MQM 
H as opposed to the boundary state approach. In fact, we will reproduce the (1, 1) ZZ 
boundary states from IMM in a very simple manner. 

5.1. sD-brane in c = 1 string theory 

In this subsection we study the closed string fields in spacetime dual to the array 
configuration in the IMM for c = 1 string. The configuration of an array of D-instantons 
separated at distance a = liwfa' (2nR = ma) is described by 



1 



X 



\ 



3 

-a 
2 



(5.1) 



The array configuration can be expressed in a gauge invariant way as 

;-l) fc / m m, ifm|fc, 
0, otherwise. 

for all k G Z, where U = e %x l R . The operator TrU k is mapped to vertex operator [21 

r(bl) 



TrtT 



(5.2) 



TrU* 



1 



,i P X + (2-\p\)4> 



(gfc/a_g-*/2) r(-| P |) 

Therefore ( |5.2| ) corresponds to the condensation of momentum modes 



p = k/R 



(5.3) 



(Vi) = (e 



(5.4) 



|) = ^j^r(O) is divergent, and since q = e 2ni / R = e 2ni / m , (q ml / 2 -q~ ml / 2 ) 
is zero. As discussed earlier (also as in pi), to make the IMM well-defined, we should add 



Formally T( 



a small imaginery part to R, say R — ► R(l — ie) so that \q\ < 1, 



qml/2 _ q-ml/2 
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At the same time, we have 

, .w 1*1 ^ (-1)' 1 



r(- 



v l-ie J {-ie)\l\ 
The factors of e cancel, and we end up with spacetime tachyon profile 



T(X, ( j>) = J2]{Vi)cos(lX)e l<l > 



i 



2 £w^ /2cos ( /x ^ (5 ' 5) 



(/! 

= Jo(2// 4 e^) + J (2// 4 e £ ^ L ) 



1=0 

4> + iX 



This result precisely agrees with the closed string fields sourced by sD-branes computed 



from ZZ boundary states |50 



In fact, we can consider a much simpler case, - a single D-instanton sitting at X = 0. 
This is described by TrU k = 1 for all k G Z in the IMM. As above, it corresponds to a 
condensation of closed string momentum modes 

(V p = k/R ) = - <T fc/2 ) = 2z sinM ^1 (5.6) 

This is nothing but the ZZ boundary state (times the Dirichlet boundary state in Euclidean 
time direction), up to a normalization factor. The agreement between IMM calculation 
and ZZ boundary state is very reminiscent to the calculation of 0. 

5.2. sD-brane in type OB string theory 

An sD-brane in type OB theory is described by an array of alternating D- and anti-D- 
instantons separated at critical distance a = -k\J2q! '. In this subsection we will set a' = 2 
for convenience. The radius is then R = 2ma/2ir — 2m, and we have q = e %a l R = e n% l ra . 
The D-instantons are located X = |a, |a, • • • , (2m — |)a, whereas the anti-D-instantons 
are located at X = |a, |a, • • • , (2m - |)a. They can be described in the followine gauge 
invariant way 

e l 2^"m, if m|/c, _ [ e~ l z™m : if m\k, 

TrU k = { TtU h = < (5.7) 

0, otherwise. I 0, otherwise. 
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From (|3.84j ), they are mapped to NS-NS and R-R vertex operators 



where p = k/R, and ,S, <S are spin fields. On RHS of ( |5.8|) we have included the corre- 
sponding leg factors. ( |5.7[ ) corresponds to condensation of closed string modes 



(V NS>P =l) = (q ml ~ q~ ml )^M(-l) 1 

, (-!)' 
~* (K|!) 2 ' 

Z+§ (-1)' 



(5.9) 



I' + fl + ^1 - |)!] 2 

where / £ Z, and we have again used the prescription R — > R(l — ze) to regularize the 
singular terms. The NS-NS and R-R scalars in spacetime are 

00 f-iV 

r(X,0)~2^i-^(2/x)'cos(/X)e^ 

= J (2v^e^) + Jo(2v^e^), (5.10) 

00 C— 1V 1 i i 

V(X, <P)~2J2 -7nWT-T^ l+ * sin ((^ + hWe^t 

1=0 1 + 2 

In particular, we find using ( |5.10| ) that (with proper normalization) 

/OO -I 
d<f>- ^e^Jo(2^e^) 
-oo ^ 



, 4,-iX , i 4,-iX , "| 1 

+ V 2/xe J (2 V 2/ue 



(5.11) 



2 



This is the conserved s-charge of the sD-branepq . 

Extending |50|1 , we can compute the closed string fields in the weak coupling region 
from the ZZ boundary state in super-Liouville theory. Our convention is b = 1, Q = 2, so 
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that c-l = 1 + 2Q 2 = 9. The disk 1-point function for the tachyon T(P) and RR scalar 
V(P) are 



^ns(P) = -iV^smhinP)^^-^)-^, 

T(-iP) 

MP) = v / 2^cosh(7rP) ^ + ^ (2^)-^. 
The closed string fields sourced by the sD-brane are (in momentum space) 

T(P,t) = -^-fe(P)Ve-(' 1+ l )fl£ (e- !,? ' + e l£f ' 
4tcE 

n>0 



1 , , COS Et 

;iPns(P)- 



4tt£ v 'sinh(a£/2) 



V(P,t) = ^^r(P) ^(-l) n e" (ri+ ^ )a£ \e~ lEt - e lEt 



n>0 



MP)- 



47r£ r ' 'cosh(a.E/2)' 
where E = \P\. Translated into Liouville coordinates, 



(5.12) 



(5.13) 



J-oo p r(-zP) 

= J (2v^e — ) + J (2^eV), 

v^t) ~ -i r *pe-»*^ v{ h+ ip \w-» (5 - i4) 

w ' J-oo p r(i-zP) 1 w 

= 2 £ ^^t( 2 ^ sinh ((^ + \)t)e^ 



which precisely agree with ( 5.10 ). 

Again, we can consider the simpler case of a single D-instanton at X = 0. This is 
described in IMM as TrU k = 1, TrU k = for all keZ.lt then follows from Q that 

r (bl) r(|p|) ^ 

which reproduce up to a normalization factor the ZZ boundary state of super-Liouville 
theory (|5.12| ). This is also a nontrivial check of the dictionary (|3.84|) . 
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6. On Black Holes in Type String Theory 

We have shown that type OA and OB string theories have the integrable structure of 
Toda lattice hierarchy. This has enabled us to prove the T-duality between type OA and 
type OB MQM at least for perturbations by purely momentum modes or winding modes. 
The matrix model deformed by modes of winding number ±1 is equivalent (in the fx — > 
limit) to the sine-Liouville theory (A/" = 2 Liouville theory in type case), which is believed 



to be dual to string theory in the Euclidean 2D black hole pl| , |3^ , p6| . The exact free energy 
of the deformed type OA and OB MQM can be computed by solving the Hirota differential 
equations. 

Let us first consider the case of type OB MQM. The winding mode perturbations, just 
as in c = 1 string, generate the Toda integrable flow. The unperturbed free energy of type 
OB MQM is perturbatively the same as that of c = 1 string, up to a redefinition a' — > 2a! 
and an overall factor of 2 (coming from doubling the fermi sea). The method used in |ZT 



to compute the free energy perturbed by A± = t±\ can be directly applied to type OB case. 
The overall factor of 2 in the unperturbed free energy is very important, since the Hirota 
differential equations are nonlinear. In the limit of large A and fx = 0, the free energy is 
now (in units a' — 1/2) 

_r_ x R + R~ 



F{\ fx = Q,R) = -(1 - R) 2 (2R - 1) — A— - - — ln(AV2i2=T) + • • • (6.1) 

where we have exhibited the genus and genus 1 terms in the expansion. The asymptotic 
radius of the Euclidean black hole in type OB string is R = a/ a' /2 = 1/2. Expanding near 
this radius, we have! 

F=-2n(R--)M-—hiM + '-- (6.2) 

_1 _ 

2 

where M oc A 1 - R . At R = 1/2, M is expected to be the mass of the black hole. In fact, 
the effective string coupling in M = 2 Liouville theory is g s X~ 2 , which is the combination 
invariant under shifting the Liouville coordinate. The mass of the black hole goes like 
M ~ l/<7sj therefore must be proportional to A 4 . This agrees with the expectation from 
the matrix model result ( |6.2| ). 

The density of states in the black hole background typically has Hagedorn growth 
behavior 

p(M) ~ M Sl e^ M (6.3) 



We are grateful to A. Adams for discussions on this point. 
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Comparing with the free energy near Tu (|6.2|), we find 



R + R- 1 5 faA . 

si + 1 = = (6.4) 

24 12 v ; 

This can be compared to the c = 1 string case, where the asymptotic radius of the 2D 
black hole is R = 3/2, and one finds s x + 1 = -(R + _R _1 )/48 = -13/288 [0. 

One can study Euclidean black holes in type OA theory in a similar way. The uncharged 
black holes should be described by a background with the lowest NS-NS winding modes 
condensed. As in fl5.8p, the perturbation of these winding modes are represented in type OA 
MQM as deformations by the operators TrO + TrO and Trfi -1 + TrO _1 . Strictly speaking, 
the coefficients of these perturbations are not the same as the time variables that generate 
the integrable flow, but essentially related to them through nonlinear bosonization maps, as 
we have seen in section 3.4. It then appears much more complicated to solve the differential 
constraints on the free energy exactly. 

However, if we are only interested in the perturbative expansion of the free energy, 
the calculation is greatly simplified. In the T-dual type OB picture, perturbatively the two 
sides of the fermi sea decouple. The deformation by NS-NS winding modes can be treated 
as independent perturbations in the two decoupled sectors. Each sector is perturbatively 
the same as c = 1 string, up to a redefinition of a' . In the end, the free energy of type 
OA theory deformed by NS-NS ±1 winding modes is simply obtained from the solution of 
||21|| for c = 1 string with a replacement a' — > a'/2, and a factor of 2 coming from the two 
sectorsB. In units with a' = 2, the answer is 

^(A, f, = 0, R) = -1(2 - R)\R -1)Aa^- *(2-R) ln ( A ^ 3T ) + ' ' ' ( 6 ' 5 ) 

Again, near the black hole asymptotic radius R = \Ja' /2 = 1, we have 

T= -2n(R- 1)M - —In MH (6.6) 

1 2 



where M oc X 2 - R ~ A 4 . The growth of the density states is of the form ( |S.3|) with the 
exponent s\ given by 

R + R' 1 1 . , 



7 Unlike the case of type OB theory, this factor of 2 trivially multiplies the answer of the 
perturbed free energy. 
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It would be interesting to reproduce the exponents (|6.4|) , (|6.7| ) from 1-loop calculations 
in M = 2 Liouville theory. One can also consider the sectors with nonzero net DO-brane 
charge q perturbed by the lowest winding modes, which should lead to charged nonextremal 
black holes. 
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